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ABSTRACT 

We report the possible optical identification of the companion to the eclipsing 
millisecond pulsar PSR J1740— 5340 in the globular cluster NGC 6397. A bright 
variable star with an anomalous red colour and optical variability which nicely 
correlates to the orbital period of the pulsar has been found close to the pulsar 
position. If confirmed, the optical light curve, reminiscent of tidal distorsions 
similar to those observed in detached and contact binaries, support the idea that 
this is the first case of a Roche lobe filling companion to a millisecond pulsar. 



Subject headings: Globular clusters: individual (NGC6397); stars: evolution 
ries: close; pulsars: individual (PSR J1740— 5340) stars: millisecond pulsar 
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1. Introduction 

The millisecond pulsar (MSP) PSR 
J1740— 5340 was discovered during a sys- 
tematic search of the globular cluster (GC) 
system for millisecond pulsars, carried out 
with the Parkes radiotelescope (D'Amico 
et al. 2001a, 2001b). The pulsar, associ- 
ated with the globular cluster NGC6397, 
is member of a binary system with a rela- 



tively wide orbit of period ~ 1.35 days, and 
it is eclipsed for about 40% of the orbital 
phase at 1.4 GHz. In a companion pa- 
per, D'Amico et al. (2001c) provide strong 
evidences that the companion star could 
be an unusual object, and give a precise 
position for the pulsar. 

We here present the results of a deep 
search of the optical companion to PSR 
J1740-5340 in the HST archive. 



1 Based on observations with the NASA/ESA 
HST, obtained at the Space Telescope Science In- 
stitute, which is operated by AURA, Inc., under 
NASA contract NAS5-26555 
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2. Observations and data analysis 

The photometric data consists of a se- 
ries of public HST exposures taken on 
March 1996 and April 1999, retrieved 
from the ESO/ST-ECF Science Archive. 
The 1999 observations consist of 116 
WFPC2 exposures with filters F555W, 
F675W, F814W and F656N (referred here 
as V555, i?675) ^814 and Ha), spanning about 
1.8 days. The 1996 observations consist 
of 55 exposures with filters F336W and 
F439W (referred here as Uzm and -B439), 
spanning 0.4 and 0.2 days respectively. 

From the accurate timing position (RAJ 
17^ 4o m 44?589; DEC J -53° 40' 40'.'9) the 
MSP in NGC6397 turns out to be approx- 
imately at 29"E and 1Q"S from the clus- 
ter center (Djorgovski & Meylan 1993). 
We used the STSDAS program metric, to 
roughly locate the MSP in the retrieved 
WFPC2 images, and it turns out to be 
within the field of view of the WF4 chip 
in both the data set. 

2.1. Astrometry 

The exact location of the MSP in the 
HST images was obtained by searching an 
astrometric solution for a wide field CCD 
image of a region around NGC 6397. We 
retrieved from the ESO Science Archive 
an image obtained on May 1999 with the 
Wide Field Imager (WFI) at the ESO 2.2m 
telescope (at European Souther Observa- 
tory, La Silla, Chile). The entire im- 
age consists of a mosaic of 8 chips (each 
with a field of view of 8' x 16') giving a 
global field of view of 33' x 34'. Only the 
chip containing the cluster center was used. 
The new astrometric Guide Star Catalog 
(GSCII) recently released and now avail- 



able on the WEB, was used to search for 
astrometric standard lying in the WFI im- 
age field of view: several hundreds astro- 
metric GSCII reference stars have been 
found, allowing an accurate absolute posi- 
tioning of the image. 

In order to derive an astrometric solu- 
tion for the WFI image we used an appro- 
priate procedure developed at the Bologna 
Observatory. The resulting rms residuals 
were of the order of ~ 0.3" both in RA 
and Dec. By using this astrometry we were 
able to accurately locate the nominal posi- 
tion of the MSP in the WFI image and in 
the HST-WFPC2 images. 

Fig. 1 shows an enlargement of a 7" x 7" 
region of the WF4 chip, centered on the 
MSP position. The 3a error circle, which 
takes into account the global error in the 
absolute positioning of the MSP, is shown. 
The global error is fully dominated by the 
uncertainty due to the astrometric proce- 
dure. One relatively bright star (A) has 
been found within the error box of the 
MSP. Two additional objects (B, C) are 
just out-side the error circle. In order 
to investigate the nature of these objects 
we performed accurate photometric analy- 
sis of the entire HST-WFPC2 data-set re- 
trieved from the Archive. 

2.2. Photometry 

The photometric reductions have been 
carried out using ROMAFOT (Buonanno 
et al. 1983), a package developed to 
perform accurate photometry in crowded 
fields and specifically optimized to han- 
dle under-sampled point spread function 
(PSF) as in the case of the HST-WF chips 
(Buonanno & Iannicola 1989). The stan- 
dard procedure described in Ferraro et 
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al. (1997) was adopted in order to de- 
rive PSF-fitting instrumental magnitudes, 
which were finally calibrated using zero- 
points listed by Holtzman et al. (1995). In 
particular we used a sample of median- 
combined images to construct reference 
Color Magnitude Diagrams (CMDs). Fig- 
ure 2 shows multiband CMDs for stars 
detected in a region of 400 x 400 pixels 
(corresponding to 40" x 40"): the location 
of the three objects are indicated. 



■Q 



E N 

V 



Fig. 1. — A portion of a median-combined 
F675W HST WFPC2 image (chip WF4) of NGC 
6397, centered near the MSP position. The re- 
gion covers about 7" x 7". Star A is the proposed 
companion to the MSP PSR J1740-5340. 

The result of the PSF fitting procedure 
for these stars have been carefully exam- 
ined by visual inspection. From this ac- 
curate photometric analysis we found that 
the bright object A lying within the er- 
ror box of the MSP has an anomalous po- 
sition in the CMD since it is located at 
the luminosity of the TO region but it 
has an anomalous red colour; the other 
two objects (B, C) are normal Main Se- 



quence stars. Individual images were in- 
stead used to check the variability of the 
objects. Objects B and C show no sig- 
nificant time variability compared to the 
measurements uncertainties. On the other 
hand, object A shows a remarkable time 
modulation (~ 0.2 — 0.3 mag) on a scale 
of several hours. This object is the vari- 
able star WF4-1 proposed by Taylor et 
al. (2001) as a BY Draconis star. 




Fig. 2. — Multiband CMDs for stars detected in 
a region of the WF4 chip (40" x 40") containing 
the MSP position. The three stars found in the 
vicinity of the MSP error box are marked with 
different symbols (a triangle for star A and squares 
for stars B and C) and labeled with their names 
in the first panel. 

2.3. Time series 

In order to check the association of the 
time variability of Star A to the pulsar bi- 
nary motion, we have carried out a period 
search analysis. The data available consist 
of four time series in the Ha, i?675, hu, 
and V555 bands taken in 1999, spanning ~ 
1.8 days, and two time series in the -B439 



3 



and {7 336 bands taken in 1996, spanning 0.2 
and 0.4 days, respectively. The periodic- 
ity search was carried out using GRATIS 
(GRaphycal Analyzer of Time Series), a 
software package developed at the Bologna 
Astronomical Observatory (see Clementini 
et al. 2000, 2001). Periods and amplitudes 
were derived for the Ha-band time series 
using GRATIS x 2 Fourier fitting routine. 
This algorithm is almost equivalent to the 
Lomb-Scargle periodogram (Scargle 1982), 
but has the advantage to have sensitivity 
also to periodicities whose light curve is not 
strictly sinusoidal, and it is more reliable 
when the data span is of length compara- 
ble to the time scale of the stellar variation 
(Faulkner 1977). As already mentioned, 
the data span of the Ha-band time series 
is ~ 1.8 days, and the searched periodic- 
ity is 1.35 days, so the x 2 fitting method is 
largely preferable. 
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Fig. 3. — Reduced \ 2 resulting from the Fourier 
fitting of the Ha data as a function of the modu- 
lation period. Small panel: light curve of the Ha 
data obtained adopting the period and the refer- 
ence epoch of the radio ephemeris and fitting the 
spectral amplitudes of the I s * and 2 nd harmonics. 

Fig. 3 shows the reduced x 2 result- 
ing from the Fourier fitting of the Ha 
data as a function of the modulation pe- 



riod. The most significant feature is indeed 
a periodicity around the predicted period 
of 1.35 days, with substantial power also 
near the 2 nd harmonic. The confidence 
level peak (99.6%) corresponds to a period 
P=1.37±0.05 days (consistent, within the 
uncertainties, with the period quoted by 
Taylor et al. (2001) for the WF4-1 vari- 
able). The quoted uncertainty corresponds 
to the period range for which the confi- 
dence level is larger than 99%, and it is 
dominated by the relatively short (~ 1.8 
days) data span available. 

Using another option of the GRATIS 
package, we have then fixed the period P 
and the reference epoch T to the radio 
ephemeris values, and have fitted the same 
Ha data for the best spectral amplitudes 
of the 1 st and 2 nd harmonics. The best- 
fit light curve, shown in the small panel 
in Fig. 3, is not exactly what we would 
expect on the basis of a simple pulsar- 
irradiation model, but as we will discuss 
in the next section, it could be understood 
in term of tidal distorsion effects occurring 
in the companion star to PSR J1740-5340. 
The time variability observed in the R, I, 
and V-bands at the same epochs follows 
a similar pattern to that observed in the 
Ha-band, and the period search analysis 
produces similar results. Fig. 4 shows 
the same 1999 Ha-band data and the t/ 336 - 
band data taken on 1996, phased using 
the accurate radio ephemeris. Remarkably, 
they show the minimum at the same or- 
bital phase, giving further evidence that 
the optical modulation is indeed associated 
to the pulsar binary motion. 
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Fig. 4. — Ha-band data and U336-band data, 
taken 3 yr a part, phased using the radio 
cphemeris. Phase 0.0 corresponds to the time of 
the ascending node. 

2.4. Is Star A the pulsar compan- 
ion? 

Can we claim that Star A (WF4-1) is 
not a BY Draconis system but is indeed 
the optical companion to the MSP? There 
is no doubt that Star A is variable, and 
that the modulation period is compatible 
to the radio orbital period. However is 
difficult to estimate the chance occurrence 
probability to find a variable star with a 
period of ~ 1.3—1.4 days in such a small 
error circle. HST observations (Taylor et 
al. 2001) have discovered several variables 
in a WFPC2/HST field of view. Also, ac- 
cording to Taylor et al. (2001), a modula- 
tion period of the order of ~day is typical 
of most BY Dra systems. 

On the other hand, in two observations 
taken three years a part, we find the min- 
ima exactly at the same phase with respect 
to the precise radio orbital period, whilst 



the light curve shape of the BY Dra sys- 
tems are expected to change, according to 
variations in the configuration of the spot- 
ted regions (Alekseev 1999) of their con- 
vective envelopes. Also, the position of 
Star A in the CMD is anomalous for a 
BY Dra or whichever other kind of binary 
system comprising two MS stars. Further 
support to the proposed association derives 
from the detection of PSR J1740-5340 in 
a Chandra pointing of NGC 6397 (Grind- 
lay et al. 2001b): its X-ray luminosity and 
color appear similar to those of the MSPs 
seen in 47 Tuc (Grindlay et al. 2001a) and 
its positional coincidence with Star A is 
consistent with the Chandra astrometric 
uncertainties. 

3. Observed properties of the com- 
panion star 

In the radio timing paper, DAmico et 
al. (2001c) demonstrate that the compan- 
ion star can not be the typical WD found 
in most binary MSPs. They propose that 
the companion can be a MS star acquired 
by exchange interaction in the cluster core 
or alternatively the same star that spun 
up the MSP and that would be still over- 
flowing its Roche lobe. Assuming that 
Star A is the pulsar companion, we here 
discuss these two hypotheses, comparing 
them with observed optical properties of 
Star A . 

In order to get some quantitative hints 
on the effective temperature T e ff and the 
radius R c of Star A , we used the recent set 
of isochrones by Silvestri et al. (1998) and 
by Vanderberg (2000). By comparing the 
CMDs in Fig. 2 with those isochrones, for 
metallicity [Fe/H]=— 2.00 and ages of t — 
12 — 14 Gyr, (compatible with the values 



5 



measured for NGC 6397), we derive R c ~ 
1.3 - 1.8 R© and T eff ~ 5500 - 5800 K for 
Star A . 

The peculiar nature of Star A can be 
unveiled inspecting the amplitude and the 
shape of its light curves (Figure 3, 4). 
There are 2 other eclipsing MSPs (PSR 
B1957+20 (Callanan, van Paradijs & Ren- 
gelink 1995) and PSR J2051-0827 (Stap- 
pers et al. 1999) ), both in the Galac- 
tic field, whose optical companion displays 
strong modulations, interpreted as due to 
the heated side of the companion entering 
in and out of view according to the orbital 
motion. Similar trend, though with a much 
smaller degree of modulation, is seen in 47 
Tuc U op t, the first identified MSP compan- 
ion in a GC (Edmonds, et al. 2001). 

The light curves of Star A are com- 
pletely different. We locate the phase 0.0 
at the ascending node of the MSP orbit; 
thus at the phase 0.75 we see the side of the 
companion facing the pulsar. In contrast 
to the other known variable MSP compan- 
ions, the light curves of Star A display 
there a minimum instead of a maximum 
(see Figure 4). Within the limits in the 
orbital period coverage of our photometry, 
the best-fit light curve of Figure 3 shows 
two maxima and two minima during each 
binary orbit: thus, tidal distorsions ap- 
pear the more natural responsible for this 
shape. They have been already invoked 
for explaining the light curves of the op- 
tical companions to black-hole candidates 
(van der Klis, et al. 1985) and NSs (Zurita, 
et al. 2000). In this scenario the maxima 
correspond to quadratures (phases 0.0 and 
0.50), when the distorted star presents the 
longest axis of its ellipsoid to the observer, 
the minima to the conjunctions. It is easy 



to recognize this trend in the insert of Fig- 
ure 3. 

Whether this is the correct interpre- 
tation, it turns out in severe constraints 
on the mass and the nature of Star A. 
The degree of ellipsoidal variations de- 
pends roughly on (Russell 1945) Am = 
k x (M M sp/M c )(R L /a) 3 F 3 sin 2 i, where M MSP 
and M c are the masses of the MSP and 
its companion, R L is the Roche lobe ra- 
dius of the companion, a is the orbital 
separation, F is the ratio between the av- 
erage radius of the star and the Roche 
lobe radius and % is the inclination of 
the orbit. The term k\ = 2.6 accounts 
for limb and gravity darkening for Ir- 
radiation from our source (Lucy 1967). 
Given the orbital parameters of PSR 
J1740-5340 (D'Amico et al. 2001c), it fol- 
lows (M MSP /M c )(R L /a) 3 ~ 0.07 for all the 
possible companions, and thus Am^ < 
0.2F 3 sin 2 i. As the observed modulation 
(Figure 4) is just ~ 0.2 mag, we argue 
that only a companion almost filling its 
Roche lobe (F ~ 1) and nearly edge-on 
(i ~ 90°) can reproduce that. Remark- 
ably, these two requirements are contem- 
porary accomplished by a companion of 
mass < 0.25 M Q , whose Roche lobe radius 
just matches the lower limit inferred for 
the observed radius R c of Star A . 

The light curve of a star affected only by 
tidally distorsion would have the minimum 
at phase 0.25 less deep than that at phase 
0.75. The reversal of this rule in the case 
of Star A can result from the overheating 
of the side facing the pulsar. In contrast 
to 47 Tuc U opt (having F ~ 0.17), and to 
the companions to PSR J2051-0827 and to 
PSR B1957+20 (for which F = 0.5 and 
F = 0.9), Star A fills up its Roche lobe 
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allowing ellipsoidal variations to dominate 
over the thermal modulation. A direct de- 
tection of both the minima would allow a 
measurement of the fraction of the imping- 
ing power from the pulsar which goes in 
heating of the surface of Star A , a very in- 
teresting value for understanding the com- 
position of the pulsar energetic flux. 

4. Discussion 

In summary, PSR J1740-5340 appears 
as the first example of a MSP orbiting 
a Roche lobe filling companion, whose 
brightness would allow unprecedented de- 
tailed investigations, for example about 
the origin of this system. 

A first hypothesis is that Star A is a 
MS star perturbed by the energetic flux 
emitted from the MSP. The so-called il- 
lumination mechanism (D'Antona 1995) 
predicts that if the heating luminosity 
L h < {l/4:)(Rja) 2 L irr (where R* is the 
star radius and L irr the MSP luminos- 
ity) is large enough, the star inflates and 
increases the effective temperature, thus 
modifying its photometric characteristics 
(Podsiadlowski 1991). The rotational en- 
ergy loss from the MSP is L{ rr ~ 1.4 x 
10 35 erg/s (D'Amico et al. 2001c). At the 
distance a ~ 6.5 R , this corresponds to 
a characteristic temperature for the heat- 
ing bath in which the star is immersed 

T h = [l/(16(J7r)L irr /a 2 ] 1/4 £ 4000 K where 
a is the constant of Stefan-Boltzmann. We 
expect that the MSP flux significantly af- 
fects the companion only if Th > T* (where 
T* is the effective temperature of an unper- 
turbed MS star) and T* < 4000 K implies 
< 0.4 M Q . As T h does not depend on R„, it 
seems energetically difficult to explain an 
increasing of ~ 40% of the effective tem- 



perature; however, only detailed simula- 
tions (Burderi, D'Antona & Burgay 2001, 
in preparation) of the system will allow to 
assess if such a low mass MS star of radius 
~ 0.2 — 0.4 R can indeed be bloated up 
to < 1.3 R and heated from ~ 4000 K 
to ~ 5500 K by the energetic flux of the 
MSP. 

Another fascinating possibility is that 
PSR J 1740 -5340 is a new-born MSP, 
the first one observed just after the end 
of the process of recycling. In this case 
Star A could have been originally a MS 
star of 1 — 2 M Q , whose evolution trig- 
gered mass transfer towards the compact 
companion, spinning it up to millisecond 
periods (Alpar, et al. 1982). Irregulari- 
ties in the mass transfer rate M c are com- 
mon in the evolution of these systems (e.g. 
Tauris & Savonije 1999): even a short 
decreasing of M c can have easily allowed 
PSR J1740-5340 (having a magnetic field 
~ 8 x 10 s G and a rotational period < 3.65 
ms) to became source of relativistic parti- 
cles and magnetodipole emission, whose 
pressure (i) first swept the environment 
of the NS, allowing coherent radio emis- 
sion to be switched on (Shvartsman 1970) 
and (ii) then kept on expelling the mat- 
ter overflowing from the Roche lobe of 
Star A (Ruderman, Shaham & Tavani 
1989). For a wide enough binary system 
(as is the case of PSR J1740-5340), once 
the radio pulsar has been switched on, any 
subsequent restoration of the original M c 
cannot quench the radioemission (Burderi 
et al. 2001). In this case we have now 
a donor star still losing matter from its 
Roche lobe at M c > 5 x 10~ n M /yr, 
(D'Amico et al. 2001c) (a high mass loss 
rate, difficult to explain in the model of a 
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bloated star). At the same time, accretion 
on the NS is inhibited due to the pres- 
sure exerted by the pulsar on the infalling 
matter. This strong interaction between 
the MSP flux and the plasma wind would 
explain also the irregularities seen in the 
radio signals from PSR J1740— 5340, some- 
times showing the presence of ionized mat- 
ter along the line of sight even when the 
pulsar is between Star A and the observer. 
The characteristic age of PSR J1740-5340 
(~ 3.5 x 10 s yr) seems indicating it is a 
young MSP, further supporting this sce- 
nario. 

If Star A will continue releasing matter 
at the present rate M c , PSR J1740-5340 
is not a candidate for becoming an isolated 
pulsar. When Star A will have shrunk well 
inside its Roche lobe, the system will prob- 
ably end up as MSP+WD (or a light non 
degenerate companion). If Star A will un- 
dergo a significant increasing of M c , the 
condition for the accretion could be re- 
established and PSR J1740-5340 would 
probably appear again as a Low Mass X- 
ray Binary or as a Soft X-ray Transient 
(Campana, et al. 1998). 
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